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It is demonstrated that the two-point correlation function of the magnetic helicity in the case of 
zero mean magnetic field has anomalous scalings for both, compressible and incompressible turbulent 
helical fluid flow. The magnetic helicity in the limit of very high electrical conductivity is conserved. 
This implies that the two-point correlation function of the conserved property does not necessarily 
have normal scaling. The reason for the anomalous scalings of the magnetic helicity correlation 
function is that the magnetic field in the equation for the two-point correlation function of the 
magnetic helicity plays a role of a pumping with anomalous scalings. It is shown also that when 
magnetic fluctuations with zero mean magnetic field are generated the magnetic helicity is very 
small even if the hydrodynamic helicity is large. Astrophysical applications of the obtained results 
are discussed. 

PACS numbers: PACS number(s): 47.65.+a, 47.27.Eq 
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I. INTRODUCTION 

Problems of intermittency and anomalous scalings for 
vector (magnetic) and scalar fields passively advected by 
a turbulent fluid flow are a subject of active research in 
the last years (see, e.g., @, § §, |, §, §, § §). The anoma- 
lous scaling means the deviation of the scaling exponents 
of the correlation function of a vector (scalar) field from 
their values obtained by the dimensional analysis. An 
interesting question is the role of the conservation laws 
in the problem of intermittency and anomalous scalings. 
For the passive scalar advected by incompressible and ho- 
mogeneous turbulent fluid flow the quantity n 2 (or T 2 ) is 
conserved (for infinitely small diffusivity or thermal con- 
ductivity), where n is the number density of particles, 
and T is the fluid temperature. Corresponding two-point 
correlation function (n(t, x)n(i, y)) has normal scaling 
(see, e.g., (§,§]). For the passive vector (magnetic field) 
the quantity h 2 is not conserved and the second moment 
(h(i,x) • h(i,y)) has anomalous scalings ||| fjj, where h 
is the magnetic field. In this case the total (magnetic 
plus hydrodynamic) energy is conserved. On the other 
hand, the magnetic helicity (&(t, x) • h(t, x)) in the limit 
of very high electrical conductivity is conserved. Here a 
is a vector potential of magnetic field, i.e., h = Vxa. 
What is a scaling for the two-point correlation function 
(a(t,x) • h(i,y)) of the magnetic helicity? 

In this Letter we show that the two-point correlation 
function of magnetic helicity has anomalous scalings for 
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both, compressible and incompressible turbulent helical 
fluid flow. For the helical fluid velocity field a(r) ^ [see 
Eq. (||) below]. The reason for the anomalous scalings 
is that the magnetic field in the equation for the two- 
point correlation function of magnetic helicity plays a 
role of a pumping with anomalous scalings. This implies 
that the two-point correlation function of the conserved 
property does not necessarily have normal scaling. This 
demonstrates a difference between passive scalar and vec- 
tor fields. We also study here an excitation of the mag- 
netic helicity by a helical turbulent fluid flow in the case 
of generation of magnetic fluctuations with zero mean 
magnetic field. Note that for a nonhelical velocity field 
[a(r) = 0, see Eq. (|J) below] the scaling exponent of the 
magnetic helicity correlation function is normal, i.e., it 
can be obtained from the dimensional arguments [w . 



II. GOVERNING EQUATIONS 

We study the evolution of magnetic fluctuations with 
zero mean magnetic field in a low-Mach-number com- 
pressible turbulent fluid flow. A mechanism of the gen- 
eration of magnetic fluctuations with a zero mean mag- 
netic field was proposed in Q and comprises stretching, 
twisting and folding of the original loop of a magnetic 
field. These non-trivial motions are three dimensional 
and result in an amplification of the magnetic field. The 
magnetic field is determined by the induction equation 
dh/dt + (u • V)h = (h • V)u - h(V • u) + f?Ah, where u is 
the fluid velocity, rj is the magnetic diffusion. We derive 
equations for the second-order correlation functions of 
the magnetic field and the magnetic helicity and we use 
a method of path integrals and modified Feynman-Kac 
formula (see, e.g., 0, H [5J ||, [n)). The equation for the 
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second-order correlation function hi 
of the magnetic field is 



{hi(t^)h 3 (t iy )) 



dhij/dt = [Lik(x)Sj a + L js (y)S l k + N ljks ]hk s + hj (I) 
(for details see [£[), where 
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and f mn = ("ru m (x)u„(y)), and S mn is the Kro- 
necker tensor and £ikm is the Levi-Civita tensor, f]ij — 
(VppSij - and U = V - Vp(ru p u)/2, and a mn = 

-(£mji{TUiV n Uj)+e n ji(TUiV m Uj))/2, &ndr/ pm = r]S pm + 
(rUpUm), and r(r) is the scale-dependent momentum re- 
laxation time, V = (u), the tensor Iij is determined by 
an external source of magnetic fluctuations, r = y x. 
We seek a solution for the second moment of the mag- 
netic field in the form 



(M*)Mx + r)) - W(r)S mn + (rW'/2)P n 



mnp 1 p 



/2 



(2) 



where P mn = 5 /r 2 , and W = dW/dr. Note 

that the current helicity correlation function fi(r) = 
[(h(x) • (Vxh(y))) + (h(y) • (Vxh(x)))]/2_ = § x /3, 
where the magnetic helicity correlation function x( r ) — 
[(a(x) • h(y)) + (h(x) ■ a(y))]/2, and S X = x" + & /r. 
The correlation function (Tu m u n ) is 

(™ m (x)ii„(x + r)) = t] T [(F + F c )8 mn + rF'P mn /2 
+ F' c r m r n /r + a{r)£ mnp r p /2] (3) 

(see |y), where ??t = UqTq/3 is the turbulent magnetic 
diffusion, uq is the characteristic velocity in the maxi- 
mum scale of turbulent motions, tq = Iq/uq, F(0) = 
1 — F c (0). The function F c (r) describes the compress- 
ible (potential) component, whereas F(r) corresponds to 
the vortical part of the turbulence. The function a(r) — 
-[(ru(x) • (Vxu(y))) + (ru(y) ■ (Vxu(x))>](6 Uq )- 1 . In 
Eqs. (g) and (J3J) the dimensionless distance r is measured 
in the units Iq. 

We use here the (^-correlated in time random pro- 
cess to describe a turbulent velocity field. Using the 
(5-correlated-in-time random process allows us to obtain 
the analytical results for the anomalous scalings of the 
two-point correlation functions of the magnetic field and 
magnetic helicity. The results remain valid also for the 
velocity field with a finite correlation time if the second- 
order correlation functions of the magnetic field and mag- 
netic helicity vary slowly in comparison to the correlation 



time of the turbulent velocity field (see, e.g., [Q). We also 
take into account the dependence of the momentum re- 
laxation time on the scale of the turbulent velocity field: 
r(k) = Tok 1 ~ p , where p is the exponent in spectrum of 
kinetic turbulent energy, k is the wave number measured 



in the units I 
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Using Eqs. (M-® we derive equations for the cor- 
relation functions of the magnetic field W(t, r) and the 
magnetic helicity x(£, r). Indeed, 

dW/dt = [W" + CW -$W)/m 

- 2(a - a(r))[i + 1 , (4) 
dn/dt = S[2(a - a(r))W + fJ,/m] , (5) 

where ao — a(r = 0) = — (tu-(Vxu)}/3 is the a-effect, / 
is an external source of magnetic fluctuations, and f /m = 
2/Rrn + 2[l-F-(rF c )']/3, and C = 4/r + m(f/m)', and 
£ = 2m(f + 2f' c )/r, and / = F + rF'/3, and f c = F c + 
rF' c /2> 7 and Rm = uqIo/i] ^ f is the magnetic Reynolds 
number, and the functions F(r) and F c (r) are determined 
below. Equation (||) and (||) are written in dimensionless 
variables: coordinates and time are measured in the units 
Iq and To, the velocity is measured in the units Uq, the 
magnetic field is measured in the units Bq. Note that in 
p"2| the system of equations which is similar to Eqs. (^) 
and (|^) was derived. However, there are mistakes in the 
equations derived in [jT2| . Since fi(r) — Sx/3, Eq. (||) 
can be rewritten as 



dx/dt = (x" + 4x'/r)/m + 6(a - a(r))W 



(6) 



Equation (@) at r = is given by (dx/dt) r=0 = 2(x" + 
4x'/ r )i , =o/R- m - This implies that in a very high electri- 
cal conductivity limit (Rm — > oo) the magnetic helic- 
ity x( r = 0) = • h(x)) is conserved. We seek a 
solution of Eqs. (Eh and (ra) in the form: W{t, r) = 
(*(r)Vro/r 2 )exp(R) and x(t, r ) = {n(r)/r 2 ) exp(R), 
where the functions Vf'(r) and n(r) are determined by 

*"/m(r) - [r + f/(r)]* = v{r)[n" - 2n/r 2 ]/^m + I , (7) 
n"/m(r) - [T + 2/{mr 2 )]n = -u(r)* , (8) 

and v(r) — 6y/m(ao — a(r)), and / = r 2 I/y/m, and 
U(r) = (C 2 + 2C' + 4^)/4m(r). We consider the case 
of small magnetic Prandtl numbers Pr m = v/r) <c f, 
which is typical for many astrophysical and geophysi- 
cal applications, where v is the kinematic viscosity. We 
choose the following model of turbulence. Incompressible 
F(r) and compressible F c (r) components in the inertial 
range of turbulence r<j < r < 1 are given by F(r) — 
(l-r^ 1 )/(l + cr), F c (r) = (l-ri- r )a/(l + a), where a 
is the degree of compressibility, q is the exponent in spec- 



trum of the function (™, 



and rd = Re 



-l/(3-p) 



P 



is the exponent in the spectrum of kinetic turbulent en- 
ergy, and Re = uqIq/v 3> 1 is the Reynolds number. Note 
that the exponent p in the spectrum of kinetic turbulent 
energy differs from that of the function (ru m u„) due to 
the scale dependence of the momentum relaxation time 
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r of the turbulent velocity field. The relation between p 
and q is q = 2p — 1 Equation (||) for a(r) =0 and 
a = was derived in |10|_. 

The solution of Eqs!~~(0) and (||) can be obtained using 
an asymptotic analysis (see, e.g., p], H m, ||). This analy- 
sis is based on the separation of scales. In particular, the 
solutions of the Schrodinger equations (Q) and (||) with a 
variable mass have different regions where the form of the 
potential U(r), mass m(r) and, therefore, eigenfunctions 
ty(r) and n(r) are different. Solutions in these different 
regions can be matched at their boundaries. The results 
obtained by this asymptotic analysis are presented below. 



III. ANOMALOUS SCALINGS 

We study a zero mode, i.e., we obtain the solutions of 
Eqs. ^> and (§) at F = 0. In this case Eqs. (@) and (|) 
are given by 



V/m{r) - U(r)^ 
k" - 2n/r 2 



I 



(9) 
(10) 



where U(r) = U(r) — 4m(ao — a(r)) 2 , and the function 
a(r) = qo(1— r 9_1 ) for < r < 1, and a(r) = for r > 1, 
and the external source of magnetic fluctuations I(r) = 
I (l-r s ) for < r < 1, andi(r) = for r > 1, and s > 0. 
The solutions of Eqs. (|j) and ( ^0| ) have three character- 
istic regions. In region I, i.e., for < r < Rm -1 "'" 1 ', 
the functions W(r) and x( r ) are given by W(r) = 
/*(1 — /SoR-mr 9-1 ), and x.( r ) = Bi + aiaoR.nir 9+1 , where 
00 =/3 m + eo(<Z~l)/(g + 2), and/, -/oRm^ 2 ^ 2 ^ 11 , 
andai = A 1 /[(q+l)(q + 4)], and f3 m = (l + qcr)/3(l + cr), 
and£ = (l + 2a)(2 + q)(q-T)/3(l + a). In region II, i.e., 
for Rm _1 " r1 ' <C r <C 1, the functions W(r) and x( r ) 
are given by 

W(r) = A 2 m 1/2 r- 3/2 coa(blnr + (fo) + W N , (11) 



x(r) = B 2 +B 3 /r 3 + a 2 a Re{ 



} , (12) 



where Wjv = -ibr 8-9 /! 2 ^^ - q/2) 2 + |6| 2 )], and 



fo 2 = 



«2 



g 2 -4 \ / 3 + a(4-q) 
4 J ^ l + qa 

(3A 2 /V2)p- 3/2 

[((l-q/2) 2 + |6| 2 )((4-q/2) 2 + |6| 2 )]V2 



Note that when q > 2 the parameter b is a com- 
plex number and Re{r~( q / 2 ~ 1+ib ^} = r -9/ 2 + 1 cos(&lnr), 
and when q < 2 the parameter b is a real number. 
For q < 2 the solution for W(r) is given by W(r) — 
m i/2 r -a/2^ 2r -|6| + yl 4r |fe|), i n region III (r > 1), the 



functions W(»") = yl 3 r (3ao cos(3ao r ) 



sin(3a r)), 



and x(r) = B^/r 3 — 6A 3 ao r_1 sin(3aor), where we take 
into account the boundary condition for the function 
x{ r )i i-e., x( r ) ~ > f° r r ~* 00 ) and a condition 
X(r) — > for o;o — > 0. Matching the functions W"(r) 



and W'(r) at the boundaries of these regions yields 
Ai ^ A 2 ^ A 3 ^ A4 ^ Iq. On the other hand, match- 
ing the functions x( r ) and x'( r ) at the boundaries of the 
regions yields B 1 ~ / aoRm ( ^ 2)/2(<?_1) for 2 < q < 3, 
and Si - I(,a„Rm^ 2+2 ' )/2(rl) for 1 < q < 2, and 
B 2 - I a for 1 < q < 3, and B 3 - 7 a Rm" (8 " 9)/2(9_1) 
for 2 < q < 3, and B 3 - I a Q Rm {q - 6+2b)/2{q - 1) for 
1 < q < 2. 

The magnetic fluctuations are excited when the mag- 
netic Reynolds number Rm > Rm cr , where the crit- 
ical magnetic Reynolds number Rm cr is found in [Q. 
For incompressible fluid a = and p = 5/3 (Kol- 
mogorov turbulence) the critical magnetic Reynolds 
number Rm' cr ' = 412, while for compressible fluid flow 
a = 0.1 the value Rm' cr ' = 740. For a larger parameter 
of compressibility the critical magnetic Reynolds num- 
ber increases sharply up to Rm (cr) ~ 10 6 (see, 0). First, 
we consider the case Rm < Rm cr , i.e., when there is 
no self-excitation of the magnetic fluctuations. In this 
case the magnetic fluctuations are sustained by an ex- 
ternal source. The first term in Eq. ( |Il| ) for the cor- 
relation function of the magnetic field W(r) in the in- 
ertial range is given by Wa ~ r~ q / 2 ~ 1 cos(61nr + tpo). 
This corresponds to the anomalous scaling of the mag- 
netic fluctuations. The normal scaling for the second 
moment of the magnetic fluctuations is given by the sec- 
ond term in Eq. (pr]): Wn ~ r 3 ~ q . The general solu- 
tion of equation for the second-order correlation function 
of magnetic field W(r) includes solutions describing the 
anomalous and normal scalings. The anomalous scaling 
Wa ~ r~ q / 2 ~ 1 cos(61nr + ipo) can be presented as the 
real part of the power-law function r e with the complex 
exponent e = —q/2 — 1 + ib(a,q). This anomalous scal- 
ing corresponds to the deviation from the condition of the 
constant flux of magnetic fluctuations over the spectrum. 
It describes the case 2 < q < 3. When 1 < q < 2 the 
anomalous exponent in a low-Mach-number compressible 
turbulent flow is real, i.e., e = —q/2 — 1 + \b(a, q)\. In 
the case of incompressible turbulent flow (a = 0) and 
1 < q < 2 this result coincides with that obtained in 
For incompressible turbulent flow and 2 < q < 3 
the anomalous scaling of magnetic field is the complex 
number e = —q/2 — 1 — ib(a — 0, q), see Jt|. 

Now we discuss solutions for the correlation func- 
tion of the magnetic helicity. Last term in Eq. ( |T^ ) 
cx i?e{r'~(' J / 2 ~ 1+lb ''} corresponds to the anomalous scal- 
ing for the magnetic helicity correlation function. The 
magnetic helicity xo — x( r — 0) in the limit of very large 
electrical conductivity is conserved. This means that the 
two-point correlation function of conserved property has 
anomalous scaling. The reason for the anomalous scal- 
ings of the magnetic helicity correlation function is that 
the magnetic field in the equation for the two-point cor- 
relation function of the magnetic helicity plays a role of 
a pumping with anomalous scalings [see Eqs. (|) and 
(|Io|)]. This demonstrates a difference between passive 
vector (magnetic field) and passive scalar. Indeed, the 
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quantity n 2 (or T 2 ) in incompressible and homogeneous 
turbulent fluid flow is conserved (for infinitely small dif- 
fusivity or thermal conductivity). Corresponding two- 
point correlation function has normal scaling. On the 
other hand, the magnetic helicity is conserved and the 
two-point correlation function of magnetic helicity has 
anomalous scaling for a turbulent helical velocity field. 
Thus we demonstrated here that two-point correlation 
function of the conserved property does not necessarily 
have normal scaling. Note that the nonlinear effects (i.e., 
self-consistent dynamics in which the back-reaction of the 
Lorentz force is considered) are important when the am- 
plitude of the magnetic field is enough large, i.e., when 
(h 2 )/47r ~ (pu 2 ). In this section we consider the case 
when there is no magnetic dynamo, i.e., the magnetic 
fluctuations are sustained by an external source /. The 
obtained results are valid when the external source / is 
not very strong, i.e., It/Ait <C (pu 2 ), and the nonlinear 
effects are not important. 

IV. DYNAMICS OF MAGNETIC HELICITY 

Now, we consider the case Rm > Rm cr , i.e., when the 
magnetic fluctuations with zero mean magnetic field are 
excited. What is the dynamics of the magnetic helicity 
in a helical turbulent fluid flow? Equations (Q) and (g) 
can be rewritten in the form (Q + V)X = TX, where X 
is the vector-column with the components X\ = ^, and 
X2 = k, and the matrix Q = Qij = 0, when i ^ j, and 
On = {d 2 /dr 2 - U)/m, and Q22 = (d 2 /dr 2 - 2/r 2 )/m, 
and the matrix V = Vij = 0, when i = j, and 
V12 = -(d 2 /dr 2 - 2/r 2 )/9m, and V21 = 1. We con- 
sider the modes which satisfy the following property: the 
change a(r) — > —a(r) in Eqs. ([?]) and (@) results in 
the change n(r) — » — re(r) and \&(r) — > ^f(r). We seek a 
solution of this equation in the form X = YlkLi x k e k + 
J y(7)E(7) dry, where the eigenfunctions e& and E(7) sat- 
isfy to the equations: Qe^ — 7fcefc, and QE(7) = 7E(7). 
Here is the vector-column with the components e±k = 
^fe, and 62k = 0, and E(7) is the vector-column with 
the components E\ = 0, and E2 = ^(7). The functions 
$ k(r) and k(r) are determined by Eqs. (Q) and (||) with 
the condition a(r) = 0. The equation for the function 
efe describes nonhelical component of the magnetic field 
correlation function and it has discrete spectrum. On 
the other hand, the equation for the function E(7) deter- 
mines helical component of the magnetic field correlation 
function and it has continuous spectrum. The continuous 
spectrum corresponds to 7 < 0, i.e., it describes the relax- 
ation of the magnetic helicity correlation function. The 
discrete spectrum of the equation for the function efc cor- 
responds to the generation of the magnetic fluctuations 
(7p > 0). The normalize conditions for the eigenfunc- 
tions efc and E(7) are given by J m(r)ete p dr = T(k)Sk p , 



and Jm(r)Et(7)E(7')dr = S(-y)5(-y - 7')- The stan- 
dard procedure used in quantum mechanics (see, e.g., 
Jl3|) yields the equations for the functions x p and 2/(7), 
i.e., (7 P - T)x p = L(y), and (r - 7)2/(7) - N(x), 
where N(x) = (1/5(7)) / v(r)m(r)k(r)x(r) dr, and 
L{y) = {l/T{p)) J v(r)^ k (r)m(r)y(r) dr, and x(r) = 
Y^h=i x k^k(r), and y(r) = J y(rf)Fih') dry' '. Now we use 
a perturbation theory (see, e.g., [p^|), i.e., we seek the 
solutions of the above equations in the form of series 
Z = Z<- 0) + eZ {Y) + £ 2 Z( 2) + . . . , where Z = x p ; y; T. 
The small parameter is e ~ Rm _ ' ,+2 " 2 ' r1 ' (see be- 
low). The perturbation theory yields r( 2fe+1 ) = 0, and 
r(°) = 7p , and r< 2 ' = -L{y^)/x p \ and y^ = 0, and 
y^ = iV(x(°))/( 7p - 7), and x ( p k+1) = 0, and x p 0) = 1, 
and Xp 2 ' is determined by equation iV(x( 2 ))(7 p — 7) = 
r( 2 )yW,etc. 

For nonhelical turbulence [a(r) = 0] the helical fi(r) 
and nonhelical W(r) parts of the magnetic field correla- 
tion function are decoupled. This implies that the mag- 
netic helicity can only relaxate from the initial value. 
In helical turbulence [a(r) 7^ 0] the magnetic helicity 
x(r) depends on W(r). This implies that the eigen- 
function e p of the discrete spectrum is modified, i.e. 
it has spiral and nonspiral components: e\ p = 
and e 2p = K p , where K p (r) = / °° h(—X, r)V p {\)/[y p + 
A] d\ + 0(e 3 ), and A = | 7 |, and * p = * p + 0(e 2 ), 
and V P (X) = (l/5(-A)) / v{r')m{r')k{-\,r)^ p {r r ) dr' . 
Using this equation we calculate the magnetic helic- 
ity x( r = 0)- The result is given by x( r = 0) ~ 
-4a T Rm" (9+2)/2(9 - 1) iy(r = 0), where W(r = 0) = 
Woexp(7 p t). For nonhelical turbulence there is only the 
relaxation of the initial magnetic helicity x(t = 0) 7^ 0. 
On the other hand, for helical turbulence the magnetic 
helicity is excited due to the growth of magnetic fluctu- 
ations. However, the magnitude of the magnetic helicity 
is very small even if the hydrodynamic helicity is large, 
i.e., x(r = 0) = -4a Q T Q W(r = 0)/Rm 13/10 , where we 
use the Kolmogorov spectrum q = 5/3. The realisabil- 
ity condition M(k) > \x{k)\k (see, e.g., Q) allows us to 
estimate the maximum possible value of the magnetic he- 
licity: Xmax ~ hWRm.- 11 ^-^, where (h 2 ) = / M(k) dk. 
Therefore, x(r = 0)/ Xm ax - (a /u )Rm- q / 2 ^ « 1. 
This means that the magnetic helicity with zero mean 
magnetic field is very small even if the hydrodynamic 
helicity is large. Only nonzero mean magnetic field can 
create large magnetic helicity (see, e.g., [p~5[| ) . 
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